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UNCLASSIEIED 


Two-Stream  Instabilities 
in  Guiding- Center  Plasmas 
for  Antihydrogen  Recombination  Schemes  ‘ 

Ronald  Stowell  and  Ronald  C.  Davidson 


Plasma  Physics  Laboratory,  Princeton  University,  Princeton,  NJ  08543 

Abstract.  Two-stream  instabilities  are  studied  analytically  in  the  guiding-center  kinetic  regime, 
which  is  chosen  in  order  that  the  results  may  be  applied  to  the  mixing  of  antiprotons  and  positrons 
preceding  antihydrogen  recombination.  The  guiding-center  kinetic  description  is  valid  for  a range 
of  parameters  which  includes  cases  in  which  magnetic  fields  are  3 — 5T;  temperatures  are  4 — 
lOK;  positron  densities  are  10^  — 10**cm“^;  and  antiproton  densities  are  lO"*  — 2x  10^ cm”^.  The 
species  occupy  long,  cylindrical  columns  coaxial  with  an  outer,  conducting,  cylindrical  wall.  Plasma 
column  radii  are  between  0.05  and  0.3  wall  radii.  A constant,  axial,  externally  generated  magnetic 
field  permeates  the  system.  Linear  stability  of  the  plasma  is  studied  as  a function  of  the  species' 
temperatures,  densities  and  column  radii;  the  mean  antiproton  velocity;  and  wavenumber.  Drifting 
Maxwellian  distribution  functions  are  considered. 


INTRODUCTION 

For  fifteen  years,  cold,  stationary  antihydrogen  has  been  sought  [1,2]  for  use  in  funda- 
mental experiments,  including  tests  of  the  equivalence  principle  of  general  relativity  [1] 
and  the  CPT  symmetry  of  the  standard  model  [3].  The  protracted  nature  of  this  effort 
is  the  result  of  technical  challenges,  such  as  deceleration  and  cooling  of  the  constituent 
species  [4,  5]  and,  nonetheless,  exceedingly  slow  recombination  rates  [6],  Only  recently 
have  cold  antiprotons  and  positrons  been  allowed  to  interact,  though  no  attempt  has  been 
made  yet  to  detect  recombined  antihydrogen  [7]. 

As  the  experimentalists  successfully  accumulate  antimatter,  collective  descriptions 
of  the  interacting  species  will  become  more  useful.  Accordingly,  plasma  phenomena 
with  which  we  are  familiar  in  other  contexts  may  appear.  In  particular,  two-stream 
instabilities,  manifest  as  the  electron-proton  (e-p)  instability  in  proton  storage  rings  [8], 
could  occur  in  antihydrogen  traps  when  the  constituent  species  are  mixed.  There  is  free 
energy  to  drive  this  instability  because  the  system  in  not  in  thermal  equilibrium  during 
mixing.  The  following  describes  preliminary  research  toward  assessing  the  possibility 
of  two-stream  instabilities  in  cold  antimatter  plasmas  and  finding  the  regimes  of  stable 
operation  of  the  traps. 
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Currently  there  are  two  distinct  collaborations  seeking  to  synthesize  cold  antihydro- 
gen at  rest  in  the  laboratory  frame.  Each  group  has  its  own  scheme  for  inducing  re- 
combination. The  ATHENA  collaboration  is  considering  a method  de,scribed  as  follows. 
Antiprotons  are  sent  into  a cloud  of  positrons  where  the  antiprotons  are  cooled  by  colli- 
sions with  the  positrons  [9].  Escaping  antiprotons  would  be  caught  electrostatically  and 
sent  back  into  the  positron  cloud.  The  ATRAP  collaboration  releases  bunches  of  cold 
antiprotons  from  an  electrostatic  well  into  a positron  cloud  [7].  As  long  as  ATHENA’s 
antiproton  densities  are  extremely  low,  a plasma  de.scription  is  not  appropriate.  However, 
this  will  change  as  ATHENA  will  ultimately  .seek  to  increase  production  of  antihydro- 
gen. A plasma  model  is  already  useful  for  ATRAP’s  parameters. 


GENERAL  THEORETICAL  MODEL 

It  may  be  worrisome  that  the  coupling  parameter  T = [e^/(47t«/3)“''^^] /kaT  is  only 
one  or  two  orders  of  magnitude  smaller  than  unity.  Even  a modest  improvement  in 
densities  or  temperatures  will  force  the  .system  into  the  regime  of  a Wigner  liquid. 
However,  for  the  parameters  presently  envisioned,  a plasma  description  of  the  antimatter 
is  valid. 

The  geometry  considered  is  shown  in  Fig.  1 . Two  cylindrical,  coaxial  plasma  columns 
of  equal  radius,  composed  of  positrons  and  antiprotons,  are  coaxial  with  an  outer, 
cylindrical,  conducting  wall.  The  antiprotons  move  through  the  cloud  of  positrons  which 
are  stationary  on  average.  Some  important  approximations  which  are  made  to  solve  the 
problem  mathematically  are  as  follows.  Bunches,  or  clouds,  of  particles  of  finite  extent 
are  treated  as  having  infinite  length.  This  is  a valid  approximation  if  we  study  the  sy.stem 
on  a scale  much  smaller  than  the  bunch  length.  For  the  present  purposes,  the  bunches 
are  also  treated  as  having  the  same  radius,  uniform  density  and  Maxwellian  unperturbed 
distribution  functions /^o.  Our  attention  is  restricted  to  azimuthally  symmetric  perturbed 
distribution  functions  f^\ . The  externally  generated  part  (j)"'  of  the  electro.static  potential 
(])  is  assumed  to  be  constant  in  .space  and  time.  Collisions  are  not  included  in  the  analysis. 
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Ill  ^ [10,  11,  12,  13],  it  can  be  seen  that,  for  relevant  temporal  and  spatial  scales,  the 
antimatter  plasmas  fall  within  the  gyrokinetic  regime,  which  is,  by  definition,  the  regime 
in  which  co/(Oc^,  (e(j)/rg+)(kiP/,p)^,  k||/kj,  and  pu/Lx  are  small  quantities  [14,  15,  16], 
Here  CO  is  any  frequency  considered,  cOcx  is  the  cyclotron  frequency  for  species  s,  (]) 
is  the  electrostatic  potential,  Ts  is  the  temperature  of  species  s,  k is  any  wavenumber 
considered,  p^^  is  the  Larmor  radius  for  species  s and  Li  is  the  perpendicular  length  of 
the  system.  The  subscripts  “1.”  and  “[|”  indicate  components  perpendicular  and  parallel 
to  the  magnetic  field  B,  respectively.  Furthermore,  the  plasmas  are  in  the  drift  kinetic 
regime,  which  is,  by  definition,  that  part  of  the  gyrokinetic  regime  in  which  k^pip  C 1. 
The  plasmas  are  in  the  guiding-center  kinetic  regime,  which  is,  by  definition,  that  part 
of  the  drift  kinetic  regime  in  which  tOpj/to^^  <C  1 for  each  species,  where  co„  and  cOp^ 
are  the  cyclotron  and  plasma  frequencies,  respectively,  for  species  s. 

From  co/cOcp  <C  1,  it  follows  that  the  E x B drift  is  much  greater  than  the  polarization 
drift.  Since  Tg+/e<\>  1,  where  ())  is  the  electrostatic  potential,  the  E x B drift  is  much 

greater  than  the  curvature  and  VB  drifts.  The  Banos  drift  [17,  18]  is  extremely  small. 

If  the  longitudinal  and  transverse  waves  are  decoupled,  we  may  further  restrict  atten- 
tion to  the  electrostatic  approximation, 

V.E=-V2(r'fq-(i)«')  =471  S ‘}s[d\fs.  (1) 

■'€{e+,p}  •' 


Here  and  (|)”'  are  the  parts  of  the  electrostatic  potential  due  to  the  plasma  and  external 
sources,  respectively,  qs  is  the  charge  of  species  s,  fs{r,\,t)  is  the  distribution  function 
for  species  s and  v is  velocity.  Since,  for  each  species,  the  collision  frequencies  are 
small  compared  to  the  plasma  frequencies  cOpj,  collisions  may  be  neglected. 

The  above  approximations  lead  to  the  collisionless  guiding-center  kinetic  equation 


dfs  [cExB 
dt  ^ B2 


nh  dvii 


= 0. 


(2) 


The  distribution  functions  fs  have  units  of  volume"*  • velocity"**.  The  quantity  nis 
is  the  mass  of  species  s.  In  cylindrical  coordinates  (r,(p,z),  with  a constant,  axial 
magnetic  field  B = Bf  and  (])''’''  independent  of  position  and  time,  Eqs.  (1)  and  (2) 
may  be  simplified.  Perturbing  the  result  about  an  equilibrium  solution  (fe+o^fpoAo) 
which  is  independent  of  tp  and  z and  applying  the  Fourier-Laplace  transform  = 

where  Im[co]  > 0,  to  the  three  first- 
order  equations  for  the  perturbed  solution  (fe+\^fp\A\)  yields  an  inhomogeneous 
initial-value  problem.  The  corresponding  homogeneous  equation  is 


0)  =0, 


(3) 


^ The  radius  of  ATRAP’s  positron  bunch  was  estimated  from  Fig.  1 of  reference  [10]  to  be  0.15mm. 
The  radius  and  length  of  ATHENA’s  antiproton  bunch  and  the  radius  of  ATRAP’s  antiproton  bunch  were 
estimated  from  ASACUSA’s  data  [1 1]  to  be  0. 1 cm,  5.0cm  and  0. 1 cm,  respectively.  All  other  data  appears 
numerically  in  the  appropriate  articles  [10,  12,  13]. 
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where  cOe  = — cEo(r)/rBo  is  the  equilibrium  ExB  rotation  frequency.  Here  the  sub- 
scripts 0 and  1 indicate  unperturbed  equilibrium  quantities  and  perturbed  quantities, 
respectively.  The  boundary  conditions  are  that  O is  zero  at  the  wall,  finite  on  the  trap’s 
axis  and  continuous  at  the  plasma  edge.  The  dispersion  relation  is  obtained  using  the 
additional  condition  found  by  integrating  Eq.  (3)  across  the  plasma  edge. 


APPLICATION  TO  COLD  PLASMAS 

By  “cold,”  it  is  meant  that  both  parallel  thermal  velocities  are  small  compared  to  the 
average  parallel  velocity  of  the  antiprotons.  In  this  special  case,  which  is  a simple  place 
to  begin, 

/iO  ==  thoH(rh  - r)  G{vj_)  8(vj  - Uj-^o) . (4) 

Here  iiso  is  the  equilibrium  density,  which  is  constant  in  space;  H is  the  Heaviside 
step  function;  rh  is  the  plasma  radiu.s,  which  is  the  same  for  both  species;  and  G(vj_) 
is  an  arbitrary  distribution  in  perpendicular  velocity  space  normalized  according  to 
2Ji/o°°dVxViG(vj.)  = 1.  If  (4)  is  used  with  f = 0,  Eq.  (3)  may  be  simplified  to 


where  kr  solves 


1 = 


, ^o(k;/'iv)/|(k;rfo)  -f  Af[(k;r/^)/o(k;rn.)  J\  (krT/,)  ^ 
^/fo(k,r„.)/o(k^rfc)-/fo(k^r/,)/o(ketv)  ^ '' Mk,rh)~ 


(5) 

(6) 


Here,  u,t  is  the  average  velocity  of  species  5;  r,,.  is  the  wall  radius;  and  /„,  7„  and  K,i  are 
Bessel  functions.  The  transcendental  equation  (6)  may  be  solved  for  kf(k;),  which,  for 
any  particular  value  of  kj,  is  quantized  with  a radial  mode  number  n.  The  solutions  for 
kr,  which  are  displayed  in  Fig.  2,  are  then  substituted  into  (5). 


kf/-fe 


4 

3 

2 

1 


FIGURE  2.  Solutions  to  (6).  For  a given  radial  mode  number  «,  the  upper  and  lower  branches  arc  for 
Oi/fK.  = 0.3  and  0.05,  respectively,  which  bound  the  range  of  values  of  of  practical  intere.st. 
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Using  Mathematica®,  the  quartic  equation  (5)  for  o)  was  solved  both  analytically  and 
numerically  with  the  same  results.  Eq.  (5)  has  four  roots.  One  root  is  unstable  for  a 
range  of  values  of  k^.  The  unstable  root,  for  physical  values  of  the  trap  parameters,  is 
shown  in  Fig.  3.  The  oscillation  frequency  Re[co]/Wpe+  and  the  growth  rate  Im[co]/o)p^+ 
are  displayed  as  functions  of  k^Up^o/(i>pe+-  As  expected  [19],  for  a cold  plasma,  the 

maximum  growth  rate  scales  as 


FIGURE  3.  The  unstable  root  of  (5)  for  n = 1 . The  left  and  right  plots  are  for  ATHENA  and  ATRAP,  re- 
spectively; ri,/r„  is  0,1  and  0.125,  respectively;  o)p^,4/o)pp  is  380  and  44,  respectively;  and  Uy;jo/(rop(;+ r/,) 
is  5.8 X lO”"*  and  8,4x  10”“*,  respectively. 


APPLICATION  TO  DRIFTING  MAXWELLIAN  PLASMAS 


Now  the  unperturbed  distribution  functions  are  allowed  to  have  nonzero  parallel  thermal 
velocities  Vf^^: 


nsoH{rb-r) 
fso  = U(vj.)  exp 


2 V? 


If  Eq.  (7)  is  substituted  into  Eq.  (3)  with  ^ = 0,  Eq.  (3)  takes  the  form 


(7) 


-I 


1 

A-nc 


+ 1 y/n y=-  exp 

V 2 


tO-k;U5^0y 

1^z5'(,vz\/2  J 


1 +erf 


. to  k^o^^o 


kjV(i2\/2  J 


— kr  + kj,  (8) 

where  Xds  = v,„/ct)pj  is  the  Debye  length  for  .species  s,  erf(z)  = {2/y/Tz)  jQdte~‘  and 
kr  solves  (6),  exactly  as  for  a cold  plasma. 

Using  Mathematica®,  the  transcendental  equation  (6)  for  kr(kj)  was  solved  numeri- 
cally, and  the  result  was  substituted  into  the  tran.scendental  equation  (8),  which  was  also 
solved  numerically.  For  ATFlENA’s  and  ATRAP’s  parameters,  it  is  found  that  there  is  a 
critical  average  velocity  for  the  antiprotons  below  which  the  system  is  stable.  In  these 
cases  and  in  all  others  examined,  it  was  found  that  the  critical  velocity  of  the  antipro- 
tons is  the  positron  thermal  velocity  times  a number  of  order  unity.  For  both  ATHENA’s 
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and  ATRAP’s  parameters,  typical  maximum  growth  rates  are  found  to  be  of  the  order  of 
Upzo/v,«  » 1. 

In  the  second  section,  terms  scaling  like  the  collision  frequencies  were  neglected  in 
favor  of  terms  scaling  like  the  plasma  frequencies.  However,  the  growth  rates  can  be 
comparable  in  size  to  which  is  about  the  same  size  as  the  largest  collision 

frequency,  the  positron-positron  collision  frequency.  This  .sugge.sts  that  collisions  play 
an  important  role  in  the  damping  or  growth  of  the  mode.  Fortunately,  to  leading  order, 
positron-positron  collisions  are  described  by  O’Neil’s  collision  operator  [20,  21].  Inclu- 
sion of  collisions  is  a logical  next  step  in  the  calculations. 
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